home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / ztrsyl.z / ztrsyl
Text File  |  1996-03-14  |  4KB  |  133 lines
  1.  
  2.  
  3.  
  4. ZZZZTTTTRRRRSSSSYYYYLLLL((((3333FFFF))))                                                          ZZZZTTTTRRRRSSSSYYYYLLLL((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      ZTRSYL - solve the complex Sylvester matrix equation
  10.  
  11. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  12.      SUBROUTINE ZTRSYL( TRANA, TRANB, ISGN, M, N, A, LDA, B, LDB, C, LDC,
  13.                         SCALE, INFO )
  14.  
  15.          CHARACTER      TRANA, TRANB
  16.  
  17.          INTEGER        INFO, ISGN, LDA, LDB, LDC, M, N
  18.  
  19.          DOUBLE         PRECISION SCALE
  20.  
  21.          COMPLEX*16     A( LDA, * ), B( LDB, * ), C( LDC, * )
  22.  
  23. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  24.      ZTRSYL solves the complex Sylvester matrix equation:
  25.  
  26.         op(A)*X + X*op(B) = scale*C or
  27.         op(A)*X - X*op(B) = scale*C,
  28.  
  29.      where op(A) = A or A**H, and A and B are both upper triangular. A is M-
  30.      by-M and B is N-by-N; the right hand side C and the solution X are M-by-
  31.      N; and scale is an output scale factor, set <= 1 to avoid overflow in X.
  32.  
  33.  
  34. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  35.      TRANA   (input) CHARACTER*1
  36.              Specifies the option op(A):
  37.              = 'N': op(A) = A    (No transpose)
  38.              = 'C': op(A) = A**H (Conjugate transpose)
  39.  
  40.      TRANB   (input) CHARACTER*1
  41.              Specifies the option op(B):
  42.              = 'N': op(B) = B    (No transpose)
  43.              = 'C': op(B) = B**H (Conjugate transpose)
  44.  
  45.      ISGN    (input) INTEGER
  46.              Specifies the sign in the equation:
  47.              = +1: solve op(A)*X + X*op(B) = scale*C
  48.              = -1: solve op(A)*X - X*op(B) = scale*C
  49.  
  50.      M       (input) INTEGER
  51.              The order of the matrix A, and the number of rows in the matrices
  52.              X and C. M >= 0.
  53.  
  54.      N       (input) INTEGER
  55.              The order of the matrix B, and the number of columns in the
  56.              matrices X and C. N >= 0.
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. ZZZZTTTTRRRRSSSSYYYYLLLL((((3333FFFF))))                                                          ZZZZTTTTRRRRSSSSYYYYLLLL((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      A       (input) COMPLEX*16 array, dimension (LDA,M)
  75.              The upper triangular matrix A.
  76.  
  77.      LDA     (input) INTEGER
  78.              The leading dimension of the array A. LDA >= max(1,M).
  79.  
  80.      B       (input) COMPLEX*16 array, dimension (LDB,N)
  81.              The upper triangular matrix B.
  82.  
  83.      LDB     (input) INTEGER
  84.              The leading dimension of the array B. LDB >= max(1,N).
  85.  
  86.      C       (input/output) COMPLEX*16 array, dimension (LDC,N)
  87.              On entry, the M-by-N right hand side matrix C.  On exit, C is
  88.              overwritten by the solution matrix X.
  89.  
  90.      LDC     (input) INTEGER
  91.              The leading dimension of the array C. LDC >= max(1,M)
  92.  
  93.      SCALE   (output) DOUBLE PRECISION
  94.              The scale factor, scale, set <= 1 to avoid overflow in X.
  95.  
  96.      INFO    (output) INTEGER
  97.              = 0: successful exit
  98.              < 0: if INFO = -i, the i-th argument had an illegal value
  99.              = 1: A and B have common or very close eigenvalues; perturbed
  100.              values were used to solve the equation (but the matrices A and B
  101.              are unchanged).
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.